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Storage of quantum information encoded into true single photons is an essential constituent of 
long-distance quantum communication based on quantum repeaters and of optical quantum infor- 
mation processing. The storage of photonic polarization qubits is, however, complicated by the 
fact that many materials are birefringent and have polarization-dependent absorption. Here we 
present and demonstrate a simple scheme that allows compensating for these polarization effects. 
The scheme is demonstrated using a solid-state quantum memory implemented with an ensemble of 
rare-earth ions doped into a biaxial yttrium orthosilicate (Y2SiOs) crystal. Heralded single photons 
generated from a filtered spontaneous parametric downconversion source are stored, and quantum 
state tomography of the retrieved polarization state reveals an average fidelity of 97.5 ±0.4%, which 
is significantly higher than what is achievable with a measure- and-prepare strategy. 

PACS numbers: 03.67.Hk, 42. 50. Ex, 42.50.Md, 42.50.Gy 



Optical quantum memories [ ] allow storage of quan- 
tum coherence and entanglement through the reversible 
mapping of quantum states of light onto matter. They 
are an essential component of long-distance quantum 
communication schemes based on quantum repeaters [2] 
and of distributed quantum networks [3]. An impor- 
tant property of quantum memories is the ability to op- 
erate with true single photons (i.e. single-photon Fock 
states). Single-photon qubits can be encoded using vari- 
ous degrees of freedoms, such as temporal modes or po- 
larization. Polarization qubits are often employed, ow- 
ing to the simplicity in performing arbitrary single-qubit 
gates and projection measurements. Yet, the storage of 
a qubit in an arbitrary state of polarization is generally 
complicated by the fact that many quantum memories 
are birefringent and have a polarization-dependent ab- 
sorption, e.g. with solid-state devices. This not only 
leads to an efficiency that depends on the state to be 
stored, but also introduces a state dependent polariza- 
tion transformation. Means to avoid this drawback were 
investigated with memory-specific approaches [4-!)] based 
on polarization-selection rules of suitably chosen atomic 
transitions. Another approach, that can in principle work 
for any type of quantum memory, is based on the embed- 
ding of two quantum memories into the arms of an in- 
terferometer [10—12], which inevitably requires long-term 
intcrferomctric stability. 

In this Letter, we present a simple and stable stor- 
age scheme for polarization qubits. We experimen- 
tally demonstrate our scheme using birefringent crystals 
doped with rare-earth ions [13]. In recent years, this 
type of quantum memory experienced rapid progress and 
key features were demonstrated, such as long storage 
times [14], high-efficiency storage [15] and large tempo- 
ral multi-mode capacity [16, 17]. Operation at the single 
photon level was demonstrated [18] using the atomic fre- 
quency comb protocol [19]. This recently lead to the stor- 



age of time-bin entangled photons generated from spon- 
taneous parametric downconversion [20, 21] and to the 
heralded creation of entanglement shared between two 
crystals [22]. Here, we experimentally show that they can 
also faithfully store polarization qubits encoded in her- 
alded single photons, despite their inherent birefringence 
and absorption anisotropy. Our scheme can in princi- 
ple work for all photon-echo based quantum memories, 
independent of the material used for the implementation. 
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FIG. 1. Scheme for compensating birefringence and absorp- 
tion anisotropy using two identical quantum memories Ma 
and Mb- (a) Light propagates along +z. The principal axes 
of the index of refraction Di and D2 of each crystal are or- 
thogonal to z, and Mb is rotated 90° around 2 with respect 
to Ma- This arrangement creates a quantum memory that 
is polarization-preserving and has a constant efficiency for all 
input polarization states, (b) Alternatively, a wave plate can 
be inserted between Ma and Mb- (c) Intensity transmitted 
along the two-memory arrangement of (a) and (b) for light 
polarized along Di or D?. The total transmission is the same 
for both components, and hence for any linear combination, 
(d) The accumulated optical phase for both polarization com- 
ponents along the two-memory arrangement. 
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Our scheme is illustrated in Fig. la using two identical 
quantum memories Ma and Mb of length L each. Both 
memories are placed along the z-axis such that their in- 
put facet is oriented to contain two of the principal axes 
of the index of refraction, which we label D± and D 2 , 
and Mb is rotated at 90° (around z) with respect to 
Ma- Light propagates along +z and hits the memories 
at normal incidence. We assume that the principal axes 
of absorption coincide with D x and D 2 . We denote a\ 
(or a 2 ) the absorption coefficient along Di (or D 2 ), and 
d\ = ol\L (or g?2 = a 2 L) the associated optical depth. 
The coincidence between the principal axes of index of 
refraction and absorption is satisfied by many types of 
quantum memories, such as rare earth-ions-doped crys- 
tals of high symmetry [23]. Hence, the two linear po- 
larizations states parallel to D\ and D 2 are eigenstates 
of each memory and form an orthonormal basis that we 
use to decompose an arbitrary polarization state. Inside 
Ma, the components will experience different absorption 
and phase shifts such that an arbitrary polarization in- 
put state undergoes a non-unitary transformation stem- 
ming from the combined effect of birefringence and ab- 
sorption anisotropy [24]. Nevertheless, the two- memory 
configuration is such that after passing through both Ma 
and Mb , both components experience the same total ab- 
sorption and global phase shift, and the arrangement is 
uniformly absorbing and polarization-preserving. Using 
a rigorous decomposition of the whole memory into in- 
finitesimal longitudinal elements, we can show that the 
memory efficiency of the forward re-emitted light can be 
the same for all polarizations, and that the polarization 
is preserved (see Appendix). The same conclusion holds 
with the configuration of Fig. lb where, instead of ro- 
tating Mb, a half- wave plate inserted between Ma and 
Mb swaps the eigenstates (here, the output polarization 
is preserved up to the swap operation). 

We note that the high fidelity of our scheme is con- 
ditioned on the interferometric stability of the distance 
between the two memories, but only on the time scale of 
the storage time. This condition is experimentally simple 
to satisfy, especially if the two memories are mounted in 
a spatial configuration yielding common mode rejection 
of vibration. In particular, this is easier to satisfy than 
the long-term stability over the ensemble of all measure- 
ments that is required for the scheme based on memories 
embedded in a Mach-Zehnder interferometer [10-12]. 

We shall now present a series of measurements made 
with a pair of 1 cm long Nd 3+ :Y2SiC>5 crystals for which 
we explicitly verified that the principal axes of the in- 
dex of refraction and of absorption coincide. We stress 
that the host biaxial crystal is highly birefringent, with 
a polarization beat length on the order of 100 um at 
883 nm [25]. First, we show that the dependence of 
the optical depth on the polarization of the light can 
be strongly attenuated. As a consequence, the efficiency 
of the quantum memory is practically constant. Second, 
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FIG. 2. Measurements of (a) optical depth and (b) mem- 
ory efficency of two Nd 3+ :Y2SiOs crystals without (•) and 
with (□) compensation scheme. Incident light is linearly po- 
larized, and the a;-axis indicates the angle of the polariza- 
tion with respect to the Di axis of the crystals, (a) Without 
compensation the optical depth varies between 2.70(1) and 
0.99(1), corresponding to propagation along the two optical 
extinction axes. With compensation the peak-to-peak varia- 
tion is reduced to 16% of the mean optical depth. Lines are 
fits to Eq. (1). (b) Efficiency of 50 ns storage measured us- 
ing laser pulses. With compensation the efficiency is almost 
independent of polarization. 



we use the quantum memory for storage and retrieval of 
heralded single photons with several polarizations, and 
obtain the fidelity of the storage process using quantum 
state tomography. 

We first measured the total optical depth d on the 
4 ^9/2 - ► 4 -F3/2 transition at 883 nm of the pair of 
Nd 3+ :Y 2 SiC>5 crystals, cooled down to 3 K, while vary- 
ing the polarization of the incident light. The optical 
depth is obtained through P out = P[ n e~ d , where P m (or 
Pont) is the optical power before (or after) the two crys- 
tals. The two crystals were initially in the configuration 
of Fig. lb, but without the wave plate between the crys- 
tals. The measured optical depth is shown in Fig. 2a. 
Without the compensation scheme, the absorption varies 
strongly with the linear polarization state, as expected. 
The overall absorption for any linear polarization can be 
calculated by decomposing the light into its components 
along the two principal axes, such that the overall optical 
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depth is given by [26], 

d = - In (e~ dl cos 2 + e~ d2 sin 2 (j>) , (1) 

where d\ (or cfo) is the optical depth for polarization 
along Di (or D2), and 9!) the angle of the linear polar- 
ization with respect to D\. A fit of the data to Eq. (1) 
results in values d x = 2.70(1) and d 2 = 0.99(1). These 
values are consistent with previous measurements on a 
single crystal [16]. We then repeated the measurement 
using the compensation scheme of Fig. lb. With compen- 
sation, the minimum and maximum values of the optical 
depth, obtained from another fit to Eq. (1), were found 
to be 1.68(1) and 1.98(1). We believe that the residual 
variation can be attributed to birefringence induced by 
the windows on the cryostat, slight misalignment of the 
crystals with respect to each other, or unequal crystal 
lengths. Additionally, it is possible that the retardancc 
and orientation of the half-wave plate deviated from the 
ideal values. 

The optical depth is important for the efficiency of the 
quantum memory. In the case of the atomic frequency 
comb, the efficiency varies as 77 oc d 2 e~ d , where d — d/F 
is the average optical depth of the comb, including a de- 
pendence on the comb finesse F [19]. This means that 
variations in the optical depth translate directly to varia- 
tions in the memory efficiency. We verified this using the 
same setup as for the measurement of the optical depth, 
but now the crystals were prepared as 120 MHz wide 
atomic frequency combs with a preprogrammed storage 
time of 50 ns (see [20, 22] for details). Using laser pulses 
of roughly Gaussian shape with 18 ns full width at half 
maximum, we measured the memory efficiency, that is, 
the ratio of the area of the retrieved pulse to that of the 
input pulse, as a function of polarization. Figure 2(b) 
shows that the efficiency follows the same trend as the 
optical depth, varying from 3% to 13% without, and from 
7% to 9.5% with compensation. 

Let us now turn to the storage and retrieval of polar- 
ization qubits. The complete setup for this purpose is 
shown in Fig. 3. The experimental cycle consists of two 
phases of 15 ms each. In the first phase the inhomoge- 
neous absorption of the two crystals is shaped into an 
atomic frequency comb by optical pumping during 11 ms 
(see [20] for details). A waiting time of 4 ms is added to 
avoid fluorescence from atoms left in the excited state. 
In the second phase the storage, retrieval and analysis 
of polarization qubits, encoded on heralded single pho- 
tons, is performed. Our photon source is described in 
detail elsewhere [20, 22]. Photon pairs are generated 
by spontaneous parametric down-conversion in a peri- 
odically poled waveguide. The signal and idler photons 
with wavelengths 883 nm and 1338 nm, respectively, are 
separated and strongly filtered to match the bandwidth 
of the quantum memories. The detection of an idler pho- 
ton heralds the presence of a signal photon, which will 
act as the polarization qubit. The signal photon then 



Heralding 
detector 




Detector 1 



FIG. 3. Experimental setup for polarization qubit tomogra- 
phy. Pairs of photons are generated by spontaneous paramet- 
ric down-conversion in a periodically poled KTP waveguide 
pumped with 4 mW of light from a continous-wave laser at 
532 nm, separated by a dichroic mirror (DM) and strongly 
filtered. The detection of a photon at 1338 nm on a su- 
perconducting nanowire single photon detector heralds the 
presence of a photon at 883 nm. The latter is in resonance 
with the quantum memories Ma and Mb and acts as a po- 
larization qubit. The initial state of the qubit is prepared 
using a fiber-based polarization controller (PC) and two wave- 
plates (A/2 and A/4). An effective half-wave plate, made of 
two quarter-wave plates between the crystals, compensates 
for their anisotropy. Finally, the polarization state of the re- 
trieved photon is analyzed using two more waveplates (A/4 
and A/2) and a polarizing beamsplitter (PBS) with a silicon 
avalanche photo diode at each output. Another polarization 
controller is used to ensure that the quantum memory is al- 
ways prepared with the same polarization, independent of the 
polarization state of the photons to be stored. 

passes through a fiber-based polarization controller, a 
half-wave plate and a quarter-wave plate to prepare the 
state of the qubit to be stored. The output of the quan- 
tum memory is sent towards a polarization analyzer con- 
sisting of a quarter-wave plate, a half-wave plate, a po- 
larizing beamsplitter and two silicon-based single photon 
detectors. The light that prepares the atomic frequency 
comb passes through the same waveplates as the single 
photons. We adjust the polarization of the preparation 
light using a fiber-based polarization controller, such that 
it always has the same polarization, independent of the 
polarization of the photons to be stored. 

To show that our compensation method allows for 
faithful storage of polarization qubits, we performed 
quantum state tomography [27] on a set of five different 
input states. The results are shown in Table I. For each 
tomographic reconstruction we performed measurements 
along the three principal axes of the Poincare sphere, 
and all detections were conditioned on the simultaneous 
detection of an idler photon. 

Based on the measured number of coincidences, we 
reconstructed the density matrix of the retrieved qubit 
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Input State 


Fidelity 


-(2) 

9 si 




99.3(6)% 


7.6(3) 


\v) 


97(1)% 


6.0(3) 


\L) 


97.7(6)% 


9.4(3) 


1+) 


95(1)% 


8.0(3) 


a\H)+p\V) 


98.7(9)% 


9.2(3) 



TABLE I. Results of the measurements performed on the 
stored and released single photons for various input states. 
The fidelity was found via tomographic state-reconstruction 
using a maximum likelihood method, and the errors esti- 
mated via Monte-Carlo simulation. Additionally, the cross- 

(2) 

correlation between signal and idler photons g si , averaged 
over measurement settings and detectors, indicates the quan- 
tum character of the process. The fifth input state is obtained 
by sending \H) through a quarter- wave plate oriented such 
that a = (1 + n/2)/2 and ,0 = 1/2. 



by maximum likelihood estimation [27]. The fidelity 
with respect to the input state is shown in Table I. We 
also estimated the uncertainty on the fidelity using a 
Monte-Carlo simulation based on the Poissonian statis- 
tics of the detection events. In all cases we find fidelities 
F > 95%, which significantly surpasses the limit of 2/3 
of any classical memory based on a measure- and-prepare 
strategy [28]. This proves that we have implemented a 
quantum memory that preserves arbitrary states of po- 
larization to a high degree. We attribute the deviation 
from unit fidelity to imperfections in the state prepara- 
tion, anisotropy compensation and analysis, caused by 
non-ideal wave plates. 

Finally, we investigated the quantum nature of the 
storage process. First, a measurement of the zero-time 
auto-correlation function of the heralded signal photon 

(2) 

(before storage) yielded <j>V < 0.06, which confirms 
the single-photon nature of the polarization qubit to be 
stored. The non-classicality of the photon retrieved from 
the memory can then be revealed with another measure- 
ment, namely the zero-time intensity cross-correlation 

(2) 

g si between signal and idler fields. Specifically, assum- 

(2) 

ing auto-correlations for signal and idler fields of 

(2) 

1 < 9x < 2, where x — s for signal and i for idler, the 
signature of non-classicality between both fields becomes 

(2) (2) 

g x si ' > 2 [29]. We measured values between g K si ' — 6.0(3) 

(2) 

and g si = 9.4(3) for all the stored polarization states, 
confirming the quantum nature of the storage and re- 
trieval process. The lowest cross-correlation value allows 
to upper bound the auto-correlation of the retrieved sig- 

(2) 

nal photon to g s ^ < 0.61(3), which is still below the 
classical threshold of 1 (the upper bound is obtained by 
assuming that the source is exactly described by a two- 
mode squeezed state [22]). We note that the relatively 

(2) 

large increase of g s ^ before and after storage is almost 
entirely due to an experimental artifact associated with 



the continuous wave operation of our source of photon 
pairs [22] , and not to a detrimental effect stemming from 
the memory itself. 

To conclude, we have experimentally demonstrated a 
scheme that allows the faithful storage of polarization 
qubits encoded into true single photons using a material 
that is birefringent and has anisotropic absorption. Wc 
note that the efficiency of a photon-echo based quantum 
memory with re-emission in the forward mode is limited 
to 54% [19, 30]. To overcome this limitation, and pos- 
sibly reach 100% efficiency, one possibility is to use the 
impedance-matched quantum memory scheme in which 
a forward-emitting quantum memory is placed between 
two-mirrors with reflectivity chosen such that all the in- 
cident light can be absorbed and re-emitted [31]. Our 
scheme thus has the potential for demonstrating a high- 
efficiency solid-state quantum memory that is compatible 
with polarization qubits. It is particularly well-suited for 
rare-earth-ions-doped crystals and greatly extends their 
range of application. For example, the storage of both 
polarization and temporal modes leads to new interesting 
possibilities, such as the quantum storage of hyperentan- 
gled photons [32]. 
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tonics, the Science and Technology Cooperation Program 
Switzerland-Russia, as well as by the European projects 
QuRep and ERC-Qore. F.B. was supported in part by 
le Fond Quebecois de la Recherche sur la Nature et les 
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We note that related results have been obtained in 
other groups [33-35]. 

Appendix 

We here present a detailed theoretical description 
of our scheme for compensating birefringence and 
anisotropic absorption in quantum memories. We rep- 
resent the polarization input state in terms of a Jones 
vector that is transformed to an output state = 
Ml'J/j), where M is a 2x2 square complex matrix. 

We use the theory derived in [ ], which treats the 
effects of birefringence (or polarization-mode dispersion, 
PMD) and polarization-dependent loss (PDL). In [24] it 
is shown that any transformation M of the input state 
induced by propagation through a memory having PMD 
and PDL effects, can be decomposed into a product of 
an unitary matrix U describing an effective PMD and an 
Hermitian positive matrix T describing an effective PDL: 

M = TU. 

In the most general case when the principal axes of 
the PDL and PMD do not coincide, T and U do not 
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commute [ ]. For rare-earth doped crystals, in most 
materials commonly used, the principal axes of the ab- 
sorption coefficient and the index of refraction do coin- 
cide due to symmetry considerations. This is the case, 
for instance, when the dopant site has an axial symmetry. 
In our crystal, Nd:Y 2 Si05, the neodymium ions are in a 
crystallographic site of low symmetry C\ and Y2SiOs is 
a biaxial crystal. The principal axes of absorption and 
birefringence could thus be oriented differently. The data 
presented in the Letter shows, however, that the principal 
axes do coincide (at least in the plane perpendicular to 
the crystallographic b-axis). We will thus consider that 
U and T are diagonal in the same basis. The principal 
axes of birefringence in Y2SiOs are denoted D\ and D 2 
in the Letter. 

The matrices T and U can then be written as 

U=(J e J>) and T=( e "J 7 ' e _° /a ), 

where d\ and d 2 are the optical depths along the princi- 
pal axes D\ and D% and <p = kAnL is the phase accu- 
mulated along axis D 2 , relative to D\, due to birefrin- 
gence An = n 2 — rii (that is the difference between the 
indices of refraction along D 2 and D\). It is now eas- 
ily shown that putting two crystals in series, where the 
second crystal is rotated 90 degrees with respect to the 
first, results in the identity transformation multiplied by 
a global polarization-independent absorption (loss) 

R(-90°)TUR(90°)TU = e -( dl+da )/ 2 e^I. 

Here R(a) is the usual 2x2 rotation matrix, a is the 
angle of rotation, and I is the identity matrix. Note 
that in terms of intensity the fraction of absorbed light 
is 1 - exp[— (di + d 2 )]. 

We have so far shown that the propagation through the 
two absorbing crystals preserves the polarization state, 
and that light is uniformly absorbed (or mapped onto) 
the crystal. For quantum memories it is of course also 
important to show that the light state can be retrieved 
efficiently and that the polarization state of the retrieved 
light is preserved as well. 

We model the quantum memory as a "sum over all 
trajectories" [19, 30, 31]. Each trajectory corresponds to 



a conversion from the input mode to the output mode 
at a particular position z in the absorbing medium. Let 
us assume that the polarization of the input mode corre- 
sponds to a polarization eigenvector of the memory, such 
that its polarization does not rotate through propaga- 
tion. Moreover, we consider the case where the input 
and output modes propagate in the same "forward" di- 
rection. In this situation, the square root of the memory 
efficiency, ^/fj, can be calculated as (see Eq. (A19) in 

[ ]) 

L 

^q = J e~ az ' 2 adze' a ^ z ^ 2 = aL £ - aL / 2 = de^ 2 , 


(A.2) 

where a is the absorption coefficient, L the memory 
length and d = aL is the optical depth. The equation 
above can be understood in the following way. For each 
trajectory there is the probability amplitude for the in- 
put photon to reach point z in the medium (e~ Q2 / 2 ), 
the probability amplitude for the photon to be absorbed 
in z and re-emitted later (adz) and the probability am- 
plitude for the output photon to reach the end of the 
medium (e~ a ( L ~ z ^ 2 ). All trajectory amplitudes can be 
added coherently in the far field, leading to the final effi- 
ciency formula for "forward" memory read-out. For the 
specific case of the atomic frequency comb protocol, the 
reasoning is the same, but the optical depth is replaced 
with d = d/F, where F is the finesse of the comb [19]. 

The well-known result of Eq. A.2 applies to a single 
polarization mode. The calculation in our case is more 
general but nevertheless straightforward. The extension 
to two polarization modes for one crystal can be done as 
follows: 

L 

M' = jT(L-z) V (L-z) ( a f a ° 2dz )nz)U(z). 



where the diagonal matrix in the middle accounts for 
the conversion efficiency for the two polarization modes. 
This equation can be simplified extensively: 



M' = 



-(L-z) ai /2 o \ / aidz o \ /g-zax/2 

e i<l> e -(L-z) a2 /2 I I o a 2 dz) I e # e -«W2 



-(L-z) ai /2 aidze -z ai /2 







e i<l> e -(L-z) a2 /2 a2dze -z a2 / 



d ie - di/2 

e^d 2 e~ d ^l 2 



The effect of both crystals can now be evaluated by sum- 



ming the amplitudes of the following two possibilities; 



6 



a photon is stored in the first crystal and propagates 
through the second crystal, or the vice versa 

R(-90°)M'R(90°)M + R(-90°)MR(90°)M' 

= e i *(di +d 2 )e- (dl+d2)/2 I. 

We conclude that the total transformation of the input 
to output mode preserves the polarization state and the 
efficiency is given by the usual formula with an effective 
optical depth d\ + d 2 , as one could have expected. The 
efficiency in this configuration would thus ultimately be 
limited to 54% for di + d 2 2 [19]. This assumes, of 
course, that there is no decoherence during the time spent 
in the memory. Decoherence can, however, be included 
by simply multiplying the equation above with the ap- 
propriate decoherence factor [30]. 

It has been shown previously that efficiencies ap- 
proaching 100% can be achieved in the so-called "back- 
ward" read out configuration [19]. Using our scheme, 
however, it is clear that one cannot in general hope to 
achieve unit efficiency in the backward direction. This is 
due to the fact that different trajectories with penetra- 
tion depth z accumulate different phases, due to the bire- 
fringence. In the special case of negligible birefringence 
(4> ~ 0), but with anisotropic absorption, one can still 
achieve unit efficiency in the backward read out config- 
uration. The calculation leading to this result is similar 
to the one presented above. 

Another way of reaching 100% efficiency is to put the 
memory in a moderate-finesse cavity, based on the con- 
cept of impedance-matching [ ]. That scheme is per- 
fectly suitable for our anisotropy compensation scheme, 
since for each cavity round trip the memory works in for- 
ward emission mode. The results presented above thus 
apply to this type of memory. There exist thus a method 
to extend the ideas presented here to unit efficiency quan- 
tum memories. 
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